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Abstract. The force on a small sphere with a plasma model dielectric function and 
in the presence of a perfectly reflecting plane is considered. The contribution of both 
the vacuum modes of the quantized electromagnetic field and of plasmon modes in the 
C*~) ■ sphere are discussed. In the case that the plasmon modes are in their ground state, 

quasi-oscillatory terms from the vacuum and plasmon parts cancel one another, leading 
a monotonic attractive force. If the plasmon modes are not in the ground state, the 
\^ , net force is quasi-oscillatory. In both cases, the sphere behaves in the same way as 

| does an atom in either its ground state or an excited state. 



1. Introduction 



As is well-known, Casimir and Polder |TJ calculated the interaction energy between an 
atom in its ground state and a perfectly reflecting plate and obtained the result 
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where z is the distance between the atom and the plate and az(i£) is the atomic 
polar izability evaluated at imaginary frequency. In the limit of large z, this potential 
takes the simple form 

which can be interpreted as due to the interaction of the atom with the quantized 
electromagnetic field. In particular, the asymptotic Casimir-Polder potential can be 
expressed as 

1 CtfO) f°° r i 
V a = --a(O) (E 2 ) = —^4 / duJ (2uj 2 z 2 -l)sm2uz+2uj zcos2luz .(3) 
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Here (E 2 ) is the shift in the mean square of the quantized electric field due to the plate, 
and the integral on ui is understood to be evaluated by use of a convergence factor. For 
example, insert a factor of e _/3 ^ in the integrand, and then take the limit (3 — > after 
evaluation of the integral. The integrand is a function which is highly oscillatory, and 
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the area under each oscillation is much larger than the final result for the integral. Thus, 
if one could slightly modify the frequency spectrum, which is the integrand in Eq. (jSJ), 
both the magnitude and sign of the result could be dramatically altered. The Casimir 
effect between parallel plates also exhibits a quasi-oscillatory spectrum [21 IH], and the 
possible effects of its modification were discussed in Ref. 

In a previous paper, Ref. the force on a small dielectric sphere in the presence 
of a perfectly reflecting plate was discussed. It was shown that the interaction of the 
sphere with the quantized electromagnetic field produces a quasi-oscillatory force whose 
magnitude can be much larger than the force associated with the asymptotic Casimir- 
Polder potential, Eq. (J2J). This result suggests a modification of the frequency spectrum 
of vacuum fluctuations. More recently, the model in Ref. [5] was extended to include 
finite reflectivity of the wall and finite temperature effects [Hj. However, Barton [2j 
has suggested that the results in Refs. [Sj and [Hj are incomplete in that they do not 
include a contribution from the quantum charge fluctuations in the sphere, and that 
this effect will cancel the quasi-oscillatory terms when the sphere is in its ground state. 
This contribution, also known as radiation reaction or the effect of plasmonic modes, 
arises from the retarded self-interaction of quantum charge fluctuations in the presence 
of the mirror. 

As will be discussed below, Barton's suggestion is correct; the quasi-oscillatory 
terms cancel for the ground state, but not for other quantum states of the plasma in 
the sphere, which is essentially a three dimensional quantum harmonic oscillator. In 
this way, the sphere behaves just as does an atom. An atom in its ground state has the 
attractive interaction found by Casimir and Polder, Eq. but an atom in an excited 
state has a quasi-oscillatory interaction potential with the plate jHj. 

In Sect. El the vacuum mode contribution found in Ref. [S] will be reviewed. The 
contribution of the charge fluctuations will be computed in Sect. EJ and the results will 
be discussed in Sect. |U 

2. Contribution of the Quantized Electromagnetic Field 

In this section, we will summarize the results of Ref. [3] concerning the force due to the 
vacuum modes, or the quantized electromagnetic field. If one has a polarizable particle 
with dynamic polarizability a(u) located a distance z from a perfectly reflecting plane, 
the force on the particle may be computed from the Maxwell stress tensor. The vacuum 
modes in the presence of the plate consist of a linear combination of an incident and a 
reflected wave. Each of these interacts with the particle and creates a radiated dipole 
field. The force on the particle arises from cross terms between (1) the reflected wave 
and the dipole field induced by the incident wave, and (2) the incident wave and the 
dipole field induced by the reflected wave. After integrating over all modes, one finds 
the net force to be 
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Here cti(u) is the real part of a(oS). The appearance of the real part of the dynamic 
polarizability is a crucial feature of this result for Fy. Equation (@J) is derived using 
a dipole approximation which requires the sphere to be small compared to both the 
distance to the plane, z, and to the characteristic wavelengths associated with ai(u). 

Now consider a sphere of radius a composed of a uniform material with dielectric 
function e(u). The complex polarizability is given by 

We will take the dielectric function to be that of the Drude model, 

*{<*) = 1 ~ , ■ v (6) 
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where u p is the plasma frequency and 7 is the damping parameter. From Eqs. (0) and 
(JHJ), we find that the real part of the polarizability is given by 

Note that although a(&) has poles only in the lower half-a; plane, its real part cti(oj) 
has poles in both the upper and lower half planes. In the context of this model, the 
sphere is required to be small in the sense that a<z and 

aw p <l. (8) 

We can evaluate the integral in Eq. (j3J) by rotating the contour of integration to 
the imaginary a>axis. However, in this process we will also acquire a contribution from 
the residue of the pole of a\{uj) at u> = Q + \ i 7, where 

6 
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Here we are primarily interested in the limit of small damping, so we will take the limit 
7 — f after the integrals are evaluated. In this limit, 

The force can be written as 

F v = J v + P Vl (11) 

where Jy is the net contribution from integrals along the imaginary axis, and Py is that 
from the pole. The explicit forms of these two contributions for 7 = are 

Jv = -4^L ^ w^w e ' { ) 



and 



a 3 u 2 



. . [2Qz(2tt 2 z 2 -3)sm(2nz) + 3(2tt 2 z 2 -l)cos(2ttz)\ .(13) 

Here the positive direction is taken to be away from the plate, so a positive value for 
the force indicates repulsion. 
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3. Contribution of the Quantum Charge Fluctuations 

In this section, we will calculate the average force on the sphere due to quantum 
fluctuations of the charge density in the sphere. Here we consider only the case where 
7 = 0, so the plasma oscillations are undamped, and have a resonate frequency at the 
pole of e(u>), that is at u = Q where Q is given by Eq. (fTUJ) . Classically, the plasma 
forms a three-dimensional harmonic oscillator with a frequency Q. When the electron gas 
is described quantum mechanically, the plasma becomes a three-dimensional quantum 
harmonic oscillator. 

First, let us consider a classical oscillating dipole in the presence of a perfectly 
reflecting plate located at z = 0. Let the instantaneous dipole moment be given by 

p(t) = p cos £lt . (14) 

The electromagnetic field radiated by this dipole is reflected by the plate and returns 
to exert a reaction force on the dipole. Let E(t) and B(f) be the reflected electric and 
magnetic fields at the location of the dipole. The resulting force on the dipole is given 
by 

F P =p j djE + pxB. (15) 

(This form can be derived directly from the Lorentz force law, and differs slightly from 
the expression for the force used in Ref . [3] , which was derived from the Maxwell stress 
tensor. However, one can show that both expressions lead to the same time-averaged 
force.) 

We are interested in the component of time-averaged force which is normal to the 
plate. A lengthy but straightforward calculation leads to the result 



PL + Pi 
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, (4ftV - 3) cos(2ftz) - 6 Slz sin(2ftz) . (16) 
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Now we need to briefly summarize the quantum mechanical description of the 
plasma in the sphere. Suppose that there are iV free electrons in the sphere, and 
that each electron is displaced by u from its equilibrium position during the plasma 
oscillations. The Lagrangian for the plasma may be written as 

£ = ijVm( u 2 - ft V) , (17) 

where m is the electron mass. Thus the system is equivalent to a single particle with 
mass Nm is a three-dimensional harmonic potential with resonant frequency uj. The 
dipole moment of the system can be expressed as 

p = iVeu, (18) 

where e is the electron charge. When the system is quantized, the Cartesian components 
of u can be expressed as 

u ' = vmm ias+a ' h (19) 
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where aj and al are the annihilation and creation operators, respectively, for quanta in 
mode j. 

We can make the link between the previous description of a classical oscillating 
dipole and a dipole undergoing quantum fluctuations by identifying (pj), the quantum 
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expectation value of p? with the time-average of the corresponding squared classical 



field. Thus, we can set 



2<P?> • 



If the oscillator is in its ground state, we have 

1 



(Px)o = (P«)( 
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(20) 
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where we used the fact that 

2 4vriVe 2 3iVe 2 , , 
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where V is the volume of the sphere. Thus in the ground state, we find that Fp = —Py. 
The force due to the averaged charge fluctuations exactly cancels the quasi-oscillatory 
term coming from the vacuum modes. The net force is then 

F = Jy + Py + Fp = Jy , (23) 

which is always attractive. 

The situation is quite different of the quantum oscillator is not in its ground state. 
In this case, we can write the net force as 

16 z 4 

+ 2Qz (4ft V - 3) sin(2ftz) 



(pV) ~ (p£>o + W ~ W>o) (80V - 3) cos(2ftz) 
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(P 2 z) - (p*)o) [(4ftV - 3) cos(2ft2) - 6 Vtz sin(2ftz) 



(24) 



Consider, for example, the case where the state of the oscillator contains one quantum 
in a particular mode. For that mode, we will have (p 2 ) = 2(p 2 )o, and F will have a 
quasi-oscillatory part which is just minus the vacuum contribution. Note that Eq. (|2*4*j) 
is equivalent to Barton's |B] Eq. (4.1). 

A case of special interest is when the oscillator is in a coherent state of each mode. 
In this case, we have 

%|V>=&#>, (25) 

where \tf)) is the quantum state of the oscillator, and the j3j are three complex numbers. 
In this case, the net force becomes 



Jy + 



2^3, 
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(Re(3 x ) 2 + (Re/3,) 2 
+ 2Qz (4ftV - 3) sin(2ftz)" 
+ 2 (Re/3 2 ) 2 [(4ftV - 3) cos(2ft2) - 6Qz sin(2ftz 



8ftV - 3) cos(2ft^) 



(26) 
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Thus if the quantum state of the plasma is a coherent state with \ j3j\ of order one, then 
there is a quasi-oscillatory force on the sphere which is of the same order of magnitude 
as the vacuum contribution, Fy. 

4. Discussion and Conclusions 

We have seen that plasmonic contributions, due to quantum charge fluctuations, give 
a significant contribution to the force on the sphere. When these modes are in their 
ground state, they exactly cancel the quasi-oscillatory terms Py arising from the vacuum 
modes. In this case, only the vacuum mode term Jy which can be written as an integral 
over imaginary frequency remains. Thus if one had ignored both the plasmon mode 
contribution and the residues of the poles of ai(a>), one would have inadvertently 
obtained the correct answer. When the plasmon modes are not in the ground state, 
the net force is quasi-oscillatory. 

In the context of atomic systems, there is an extensive literature on the relative 
effects of vacuum modes and quantum charge fluctuations, often called radiation reaction 
or self-reaction. See, for example, Ref. jU] and ^U]. However, many of the treatments for 
atomic systems do not make it clear whether the separation between vacuum and self- 
reaction effects is gauge invariant. In the present model, it is obvious that the two effects 
are separately gauge invariant, as all of the calculations involve field strengths rather 
than potentials. When calculations are performed by doing second order perturbation 
theory on a Hamiltonian, as in Refs. pQ and jH], there is no need to consider the effects 
separately, and the cancellation is automatic. 

In the case of the Casimir effect between planar geometries, the role of plasmonic 
contributions at short distances has long been recognized |TT| Il2j . More recently, 
Intravaia and Lambrecht ^3] have argued that they are significant at all separations. 

It seems to be desirable to understand better the physical reason for the cancellation 
between vacuum and plasmonic effects. It is clear that one way to upset the cancellation 
is to excite the plasmon modes. However, it is less clear if this is the only way in more 
general geometries. 
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